A variational coupled-cluster study of magnon-density-wave excitations in quantum 
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We extend recently proposed variational coupled-cluster method to describe excitation states of 
quantum antiferromagnetic bipartite lattices. We reproduce the spin-wave excitations (i.e., magnons 
with spin ±1). In addition, we obtain a new, spin-zero excitation (magnon-density waves) which has 
been missing in all existing spin-wave theories. Within our approximation, this magnon-density- 
wave excitation has a nonzero energy gap in a cubic lattice and is gapless in a square lattice, similar 
to those of charge-density-wave excitations (plasmons) in quantum electron gases. 
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From many-body physics view, energy spectra of low- 
lying excitation states of a quantum many-body interact- 
ing system are mainly determined by the correlations in- 
cluded in its ground state. Feynman's description of the 
phonon-roton excitations in a Helium-4 quantum liquid is 
a well-known example 0, Q ■ For the ground state prop- 
erties, the method of correlated basis functions (CBF) 
has produced satisfactory results. Combining with Feyn- 
man's excitation theory, we therefore have an analytical, 
systematic method capable of describing most states of 
the Helium-4 quantum liquid at low temperature Q- 

The coupled-cluster method (CCM) HUIl] is another 
successful microscopic many-body theory which has pro- 
duced many most competitive results for the ground state 
properties of many atoms, molecules, electron gas, and 
quantum spin systems |6|]. In the traditional CCM, how- 
ever, the bra and ket states are not hermitian to one an- 
other • We recently extended the CCM to a variational 
formalism 0, 0] , in which bra and ket states are hermi- 
tian to one another. In our analysis, a close relation with 
the CBF method was found and exploited by using dia- 
grammatic techniques. In application to quantum anti- 
ferromagnets, our calculations in a simple approximation 
reproduced the ground state properties of the spin-wave 
theory (SWT) [1J|. Improvements over SWT by higher- 
order calculations were also obtained ; 9j . Here we extend 
this variational method to describe excitation states. In- 
spired by the close relation to the CBF method, we are 
able to obtain the energy spectra of quasiparticle-dcnsity- 
wave excitations by following Feynman, as well as the 
usual quasiparticle excitation spectra as in the traditional 
CCM. For the antiferromagnetic models, these quasipar- 
ticle excitations correspond to the well known spin-wave 
excitations (magnons) with spin ±1; the new, spin-zero 
excitations are identified as magnon-density-wave exci- 
tations whose energy spectra share some characteristics 
with that of charge-density- wave excitations (plasmons) 
in quantum electron gases [llj . We like to emphasize that 
these spin-zero excitations have been missing in all spin- 
wave theories, includi ng t he more recently modified spin- 
wave theories HI HQ Ma gnon-density fluctuations 



in Heisenberg ferromagnets at low temperature were in- 
vestigated in the sixties by calculating the longitudinal 
susceptibility [lil Hfil ] . Calculations for the longitudinal 
spin fluctuations due to multi-spin-wave excitations in 
antiferromagnets were also performed by a bosonization 
scheme |l7l| and by finite-size (classical) dynamic simula- 
tions [13 . Both our analysis and our results are different 
to these calculations and more discussion will be given in 
conclusion. 

The Heisenberg Hamiltonian is given by 
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where the index I runs over all N bipartite lattice sites, n 
runs over all z nearest-neighbor sites, s z and s are spin 
operators. As shown in Refs. 8 and 9, we use Coester 
representation for the ground-state \^> g ) of Eq. (1), 
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with its hermitian conjugate (^> g \ — (Q\e s , S = ^ 7 FiCj 
as the bra ground-state. In Eq. (2), the model state |$) 
is given by the Neel state with alternating spin up and 
down sublattices, and C\ with nominal index / is given 
by the spin-flip operators, 



N/2 

E 

fc=l i 



E 

— ,h- 



fii ■ 



3 k 



(28)* 



(3) 



with s as spin quantum number. The bra state operators 
are given by the corresponding hermitian conjugate of 
Eq. (3), using notation Fx — fi 1 ...j 1 ... for the bra-state 
coefficients. As before, we have used index i exclusively 
for the spin-up sublattice of the Neel state and index j for 
corresponding the spin-down sublattice. The coefficients 
{Fi,Fi} are then determined by the usual variational 
equations as 
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The important bare distribution functions, gj = (Cj) 
and gj = (Cj), can be expressed in a self-consistency 
equation as 



9i = G{gj,Fj), gi = G(gj, Fj), 
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where G is a function containing up to linear terms in gj 
(or gj) and finite order terms in Fj (or Fj). In Refs. 8 
and 9, as a demonstration, we considered a truncation 
approximation in which the correlation operators S and 
S of Eqs. (2-3) retain only the two-spin-flip operators as 
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The spontaneous magnetization (order-parameter) in 
this two-spin-flip approximation is obtained by calculat- 
ing the one-body density function as 



<*<} =s~Pi 
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where Pi = p due to translational invariance of the lattice 
system. As demonstrated in Ref. 9, contributions to one- 
body density function of Eq. (7) can be represented by 
diagrams. The results of SWT are reproduced by resum- 
ing all diagrams without any so-called Pauli lines repre- 
senting Pauli exclusion principle manifested by (s ± ) 2 = 
for s = 1/2. For example, the one-body bare distribu- 
tion function g^j = (s^s^)/(2s) in this approximation is 
given by, after the sublattice Fourier transformation, 
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where f q is the Fourier component of correlation coef- 
ficient fij with q restricted in the magnetic zone etc. 
Variational Eq. (4) reproduces the SWT result for this 
coefficient as 



fq fq 



f 

lq 



1 - (7 9 ) 2 - 1 



7« 



Finally, the two-body distribution functions gij,vy 
(s^ s+ sji s~j,) I {2s) 2 is approximated by 

9ij,i'j' ~ 9ij9i'j' i 9ij'9i'j ■ 



(9) 
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To go beyond SWT, we need to consider Pauli principle 
as mentioned earlier by including diagrams with the so- 
called Pauli lines. Details were described in Ref. 9. In 
the followings, we will use approximations Eq. (6-10) to 
discuss excitation states. 

Order-parameter of Eq. (7) can also be calculated 
through two-body functions as 



(11) 



where s z a — ( — l) z sf /AT is the staggered spin operator 
[l3j . Eq. (11) is in fact the sum rule for the two-body 
distribution function as in CBF 2]. This can seen by 
introducing (staggered) magnon-density operator hi as 



2hi = 2s — S4 
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where as before summation over n is over all z nearest 
neighbors. Hence, the sum rule for the one-body function 
is simply j| Y^ii^i) = P- The two-body Eq. (11) can now 
be written in the following familiar sum rule equation as 



2 N/2 

jr^2(hihi>) = ppi = p 2 , 

i' = l 



(13) 



where in the last equation, translational invariant prop- 
erty pi — p has been used. In the approximation of 
Eq. (6-10), we find that this sum rule is obeyed in both 
cubic and square lattices in the limit N — > oo. In partic- 
ular, we find that (A J^i'i^i^i') ~ P 2 ) 00 V-^ m a cubic 
lattice and oc (\nN)/N in a square lattice [l^. These 
asymptotic properties are important in the correspond- 
ing excitation states as discussed later. However, Eq. (13) 
is violated in the one-dimensional model, showing the de- 
ficiency of the two-spin-flip approximation of Eq. (6) for 
the one-dimensional model. We therefore leave more de- 
tailed discussion and a possible cure elsewhere 0] and 
focus on the cubic and square lattices in the followings, 
using approximations of Eqs. (6-10). 

For the quasiparticle excitations, briefly, we follow Em- 
rich in the traditional CCM |20j, |2l| and write excitation 
ket-state |* e ) involving only Cj operators as 



(14) 



but, unlike the traditional CCM, our bra excitation 
state is the corresponding hcrmitian conjugate (^ e \ = 
(& g \X — ($\e s X. Choosing a single spin-flip operator 
C\ = s~ , we have X w Y^i x i s 7 with coefficient cho- 
sen as Xi — Xi(q) — y^e lq ' ri to define a linear mo- 
mentum q. The state of Eq. (14) has therefore spin 
s totai = — !• The energy difference between the above ex- 
citation state and the variational ground state of Eqs. (2- 

4 ). (*al*ff*l* fl >/(*e|* e ) - ( H ), can be derived 

as [l9|, to the order of (2s), 
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where p q — f q g q and we have used Eqs. (6-10). This 
agrees with SWT 10] in this order. Spectrum of Eq. (15) 
is gapless in all dimensions because e q oc q as q — > 0. Sim- 
ilar calculations using spin-flip operators C\ — sj of the 
j-sublattice in Eq. (14) will produce the same spectrum 
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as Eq. (15) except that the corresponding excitation state 
has spin s 



total 



hi. These excitations are often referred 
to as magnons [lo|. 

We next consider the magnon-density-wave excita- 
tions. Quasiparticles such as magnons in general interact 
with one another, thus producing quasiparticle density 
fluctuations. The interaction potential for magnons is de- 
scribed by the term sfsf +n in the Hamiltonian of Eq. (1). 
Following Feynman as in CBF Q, Q , we write our gen- 
eral quasiparticle-density-wave excitation state as, using 
magnon-density operator as defined in Eq. (12), 

|*°} = X°|* S >, X? = 5><(?)rii, q>0 (16) 



and its hermitian counterparts for the bra state, (^°| = 
(V g \X°. The coefficient x t (q) = y^e lq ri as before. The 
condition q > in Eq. (16) ensures the orthogonality 
between this excited state with the ground state. We 
notice that the density operator hi in Eq. (16) is a her- 
mitian operator. This property can be used to derive 
a double commutation formula for the energy difference 
between the above excitation state and the variational 
ground states of Eqs. (2-4) as, 



„_ (V g \X°HX^ g ) 
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where N(q) = {[X°, [H,X$\])/2, S°(q) = (X°X° q ) is the 
structure function, and notation (• ■ •} is the ground state 
expectation value as before. Details will be given else- 
where [lflj . The double commutator N(q) and structure 
function S°(q) can be straightforwardly calculated as, us- 
ing approximations of Eqs. (6-10), 
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and 



S°(q) = j(l+7 9 2 )p+^[(l+7 g 2 )p 9 'P g - 9 '+27 9 5 9 '3 9 -g'], 
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where p q = f q g q as before, f q and g q are as given by 
Eqs. (8) and (9). Substituting Eqs. (18) and (19) into 
Eq. (17), we can then calculate energy spectrum e q nu- 
merically. We notice that Eq. (19) is closely related to 
the the sum rule Eq. (13) which corresponds to S°(q) 
at q = 0. It is not difficult to see from Eq. (18), N(q) 
has a nonzero, finite value for all values of q. Any spe- 
cial feature such as gapless in the spectrum e q there- 
fore comes from the the structure function of Eq. (19), 
and hence is determined by the asymptotic behaviors 
of the sum rule Eq. (13) mentioned earlier. For a cu- 
bic lattice, we find that the spectrum e q has a nonzero 
gap everywhere. The minimum gap is e q w 0.99sz at 



q — > 0. This gap is about the same as the largest 
magnon energy, e q = sz at q = (tt/2, tt/2, tt/2) from 
Eq. (15). At q = (tt/2, tt/2, tt/2), we have the largest en- 
ergy e q w 2.92sz. This is nearly three magnons' energy 
at this value of q. At q = (tt, 0, 0), we obtain e q ~ 2.56sz. 

For the square lattice the structure function S°(q) of 
Eq. (19) has a logarithmic behavior Inq as q — » 0. This 
is not surprising as discussed earlier in the sum rule 
Eq. (13), where occurs the asymptotic behavior (In N)/N 
as N — > oo. For small values of q, N(q) approaches to a 
finite value, N(q) » 0.275sz as q — > 0. The correspond- 
ing energy spectrum of Eq. (17) is therefore gapless as 
q — > 0. We plot the spectrum for the values of q be- 
tween (0.017T, 0) to (tt, 0) in Fig. 1, together with the 
corresponding magnon energies for comparison. As can 
be seen from Fig. 1, magnon-density-wave energy is al- 
ways larger than the corresponding magnon energy. At 
small values of q (q < 0.057r), we find a good approxima- 
tion by numerical calculations for the structure function, 
S°(q) w 0.24-0.16 In q. The energy spectrum of Eq. (17) 
can therefore be approximated by 
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for a square lattice. This spectrum is very "hard" when 
comparing with the magnon's soft mode e q oc q at small 
q. For example, we consider a system with lattice size of 
N = 10 10 , the smallest value for q is about q « 10~ 10 7r 
and we have energy e° « 0.07sz. Comparing this value 
with the corresponding magnon energy e q w 10~ 10 sz, we 
conclude that the energy spectrum of Eq. (17) is "nearly 
gapped" in a square lattice. We also notice that the 
largest energy in a square lattice e q w 2.79sz at q = 
(tt, 0), not at q = (tt/2, tt/2) as the case in a cubic 
lattice. At q = (tt/2, tt/2), we obtain e q w 2.62sz for the 
square lattice. 

In summary, we have extended recently proposed vari- 
ational CCM to describe excited states and have applied 
to the well-studied Heisenberg antiferromagnets. We 
have obtained both the well-known quasiparticle exci- 
tations (magnons) by using operators a and the new 
quasiparticle-density-wave excitations (magnon-density 
waves) by using operators s z . The special features and 
numerical values discussed above for the cubic and square 
lattices clearly indicate that this magnon-density-wave 
excitation can not be understood as a simple sum of 
a pair of spin ±1 magnon excitations; they are differ- 
ent to the the multi-magnon excitations as discussed in 
Refs. 17 and 18. In fact, magnon-density-wave excita- 
tions share some characteristics with charge-density- wave 
excitations (plasmons) in quantum electron gases, with 
a large energy gap in three dimension and the gapless 
spectrum in two dimension; it is also interesting to know 
that the plasmon spectrum in two-dimensional system is 
also much "harder" than the corresponding quasiparticle 
excitations ^lj. We also notice that recently modified 
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spin-wave theories (SWT) were applied to finite systems 
with results in reasonable agreements with exact finite- 
size calculations 0,0,0]. As pointed out in Ref. 13, 
however, a major deficiency in this modified SWT is the 
missing spin-zero excitations as the low-lying excitations 
for a finite lattice Heisenberg model is always triplet with 
spin equal to 0, ±1. We believe our magnon-density-wave 
excitation as discussed here corresponds to the missing 
branch; the energy gap in the cubic lattice and the nearly 
gapped spectrum in the square lattice of Eq. (17) reflect 
the nature of long-ranged Neel order in the ground states 
of infinite systems. In a recent experiment on the an- 
tiferromagnet MnF2 longitudinal spin fluctuations were 
observed and explained as two magnon excitations [22j |. 
It will be interesting to see further experiments on these 
spin- zero fluctuations in both three and two dimensional 
systems at very low temperature and high energy for pos- 
sible observation of magnon-density-wave excitations as 
described here. 
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Fig. 1 Excitation energy spectra for the values of 
q = (0.017T, 0) to (-7T, 0) in a square lattice. The higher 
branch is for the magnon-density-wave excitation of Eq. (17) 
and the lower one is for the magnon excitation of Eq. (15). 



